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Abstract 

For an rit transmit, jir receive antenna system (rit x system), a full-rate space time block code 
(STBC) transmits rimin = min{nt,nr) complex symbols per channel use and in general, has an ML- 
decoding complexity of the order of A/"'"""" (considering square designs), where M is the constellation 
size. In this paper, a scheme to obtain a full-rate STBC for 2° transmit antennas and any n^, with reduced 
ML-decoding complexity of the order of A/"''^"'"'"^^), is presented. The weight matrices of the proposed 
STBC are obtained from the unitary matrix representations of a Clifford Algebra. For any value of rir, 
the proposed design offers a reduction from the full ML-decoding complexity by a factor of i\/ * . The 
well known Silver code for 2 transmit antennas is a special case of the proposed scheme. Further, it is 
shown that the codes constructed using the scheme have higher ergodic capacity than the well known 
punctured Perfect codes for < rit- Simulation results of the symbol error rates are shown for 8x2 
systems, where the comparison of the proposed code is with the punctured Perfect code for 8 transmit 
antennas. The proposed code matches the punctured perfect code in error performance, while having 
reduced ML-decoding complexity and higher ergodic capacity. 



L Introduction and Background 

Complex orthogonal designs (CODs) [1], [2], although they provide linear Maximum Likelihood (ML) 
decoding, do not offer a high rate of transmission. A full-rate code for an rit x MIMO system transmits 
min{nt, Ur) complex symbols per channel use. Among the CODs, only the Alamouti code for 2 transmit 
antennas is full-rate for a 2 x 1 MIMO system. A full-rate STBC can efficiently utilize all the degrees 
of freedom the channel provides. An increase in the rate also results in an increase in the ML-decoding 
complexity. The Golden code [3] for 2 transmit antennas is an example of a full-rate STBC for any number 
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of receive antennas. Until recently, the ML-decoding complexity of the Golden code was reported to be 
of the order of M^, where M is the size of the signal constellation. However, it was shown in [4], 
[5] that the Golden code has a decoding complexity of the order of NP^/M only. A lot of attention 
is being given to reducing the ML-complexity of full-rate codes. Current research focuses on obtaining 
high rate codes with reduced ML-decoding complexity (refer to Sec. |ll] for a formal definition), since 
high rate codes are essential to exploit the available degrees of freedom of the MIMO channel. For 2 
transmit antennas, the Silver code [6], [7], is a full-rate code with full-diversity and an ML-decoding 
complexity of order for square QAM. For 4 transmit antennas, Biglieri et. al. proposed a rate-2 
STBC which has an ML-decoding complexity of M^a/M for square QAM without full-diversity [8]. It 
was, however, shown that there was no significant reduction in error performance at low to medium SNR 
when compared with the previously best known code - the DjABBA code [6]. This code was obtained by 
multiplexing Quasi-orthogonal designs(QOD) for 4 transmit antennas [9]. In [5], a new full-rate STBC 
for 4 X 2 system with full diversity and an ML-decoding complexity of M^\/M was proposed. This 
code was obtained by multiplexing the coordinate interleaved orthogonal designs (CIODs) for 4 transmit 
antennas [10]. These results show that codes obtained by multiplexing low complexity STBCs can result 
in high rate STBCs with reduced ML-decoding complexity and without any significant degradation in 
the error performance when compared with the best existing STBCs. Such an approach has also been 
adopted in [11] to obtain high rate codes from multiplexed orthogonal designs. 

In general, it is not known how one can design full-rate STBCs for arbitrary number of transmit and 
receive antennas with reduced ML-decoding complexity. Such a design has been presented for = 4 in 
[12]. It is known how to design information lossless codes [13] for the case where n,. > n^. However, 
it is not known how to design information lossless codes when Ur < rit. In this paper, we design codes 
which have higher ergodic capacity at high signal to noise ratio (SNR) than the best existing codes (the 
Perfect codes [14]) for < nt. The resulting codes also have lower ML-decoding complexity than the 
comparable punctured Perfect codes. The contributions of the paper are: 

1) We analyze the ergodic capacity of MIMO channels with space time codes when < nt. We 
relate the entries of the R-matrix of the equivalent channel matrix to ergodic capacity at high SNR. 

2) We give a scheme to obtain rate-1, 4-group decodable codes (refer Section|lI]for a formal definition 
of multigroup decodable codes) for nt = 2^ through algebraic methods. The speciality of the 
obtained design is that it is amenable for extension to higher number of receive antennas, resulting 
in full-rate, reduced ML-decoding complexity codes for any number of receive antennas, unlike 
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the previous constructions [15], [16], [17] of rate-1, 4-group decodable codes. 
3) We propose a scheme to obtain full-rate, reduced ML-decoding complexity codes for 2" transmit 
antennas and any number of receive antennas. These codes are also shown to have higher ergodic 
capacity than the comparable punctured Perfect codes for the case n,. < nt, and lower ML-decoding 
complexity as well. In terms of error performance, the proposed codes have more or less the same 
performance as the corresponding punctured Perfect codes. This is shown through simulation results 
for the 8 X 2 MIMO system. 
The paper is organized as follows. In Section |II1 we present the system model and the relevant 
definitions. The ergodic capacity analysis is presented in Section JII] and the method to construct Rate-1, 
4-group decodable codes is proposed in Section JV] The scheme to extend the code to obtain full-rate 
STBCs for higher number of receive antennas is presented in Section jV] Simulation results are discussed 
in Section |Vl] and the concluding remarks are made in Section IVIII 

Notations: Throughout, bold, lowercase letters are used to denote vectors and bold, uppercase letters 
are used to denote matrices. Let X be a complex matrix. Then, and X^ denote the Hermitian and the 
transpose of X, respectively and j represents The (i, j)*'^ entry of X is denoted by X(i, j) and tr(X) 

denotes the trace of X. The set of all real and complex numbers are denoted by M and C, respectively. 
The real and the imaginary part of a complex number x are denoted by xj and xq, respectively. ||X|| 
denotes the Frobenius norm of X, and and Or denote the T x T identity matrix and the null matrix, 
respectively. The Kronecker product is denoted by ^. For a complex random variable X, £[X] denotes 
the mean of X. The inner product of two vectors x and y is denoted by (x,y). 
For a complex variable x, the (.) operator acting on x is defined as 



~ A 
X = 



XI -XQ 
XQ XI 



The (.) can similarly be applied to any matrix X G Qnxm j-gpiacing each entry Xij by Xij, i 
1, 2, • • • , n, j = 1, 2, • • • , m , resulting in a matrix denoted by X e M2nx2m 
Given a complex vector x = [xi,X2, - ■ ■ , Xn]'^, x is defined as 



X = [xil,XiQ, ■■■ , Xnl, XnQ]"^. 
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II. System Model 

We consider Rayleigh block fading MIMO channel with full channel state information (CSI) at the 
receiver but not at the transmitter. For nt x MIMO transmission, we have 

fSNR 

Y=J HS + N, (1) 

V nt 

where S G C"'^^ is the codeword matrix whose average energy is given by <5(||S|p) = ?ifT, transmitted 
over T channel uses, N G C"-^^^ is a complex white Gaussian noise matrix with i.i.d entries ~ Ac (0, 1) 
and H G xnt jj^g channel matrix with the entries assumed to be i.i.d circularly symmetric Gaussian 
random variables ~ A/c (0, 1). Y G c^^xT ^j^g received matrix and SNR is the signal to noise ratio 
at each receive antenna. 

Definition 1: (Code rate) Code rate is the average number of independent information symbols trans- 
mitted per channel use. If there are k independent complex information symbols (or 2k real information 
symbols) in the codeword which are transmitted over T channel uses, then, the code rate is k/T complex 
symbols per channel use {2k /T real symbols per channel use). 

Definition 2: (Full-rate STBCs) For an nt x MIMO system, if the code rate is min{nt,nr) 
complex symbols per channel use, then the STBC is said to he full-rate. 

Assuming ML-decoding, the ML-decoding metric that is to be minimized over all possible values of 
codewords S is given by 

M(S) = ||Y-./^HSf (2) 

V nt 

Definition 3: (ML-Decoding complexity) The ML decoding complexity is measured in terms of the 
maximum number of symbols that need to be jointly decoded in minimizing the ML decoding metric. 
For example, if the codeword transmits k independent symbols of which a maximum of p symbols need 
to be jointly decoded, the ML-decoding complexity is of the order of AP, where M is the size of the 
signal constellation. If the code has an ML-decoding complexity of order less than M'^, the code is said 
to admit reduced ML-decoding. 

Definition 4: (Generator matrix) For any STBC that encodes 2k real symbols (or k complex infor- 
mation symbols), the generator matrix G is defined by the following equation [8]. 

vec (S) = Gs, 

where S is the codeword matrix, s = [si, S2, • ■ ■ > •S2fc]^ is the real information symbol vector. 

A codeword matrix of an STBC can be expressed in terms of weight matrices (linear dispersion 
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matrices) [18] as 

2k 



S — ^ ^ -SjAj. 



i=l 

Here, Aj,i = 1,2,- •• ,2k are the complex weight matrices for the STBC and should form a linearly 
independent set over M. It follows that 

G = [uec(Ai) uec(A2) ••• fec(A2A;)]. 

Definition 5: (Multigroup decodable STBCs) An STBC is said to be 51-group decodable [17] if its 
weight matrices can be separated into g groups Qi, Q2, ■ ■ ■ , 9g such that 

AiAf + AjAf = On, , Aiegi,Aj €gp,l^p. (3) 
Equation ([T|| can be rewritten as 



SNR^^ -~r, 

vec(Y) = \ He„s + vec{H), 

V rit 

where Hgg G ]^2n,rx2n„,„T given by 

Heg = (It ® H) G, 
with G G ]^2ntTx2n„,„T (j^g generator matrix as in Def. HI 

III. Relationship between weight matrices and ergodic capacity 

It has been shown that if the generator matrix is unitary, the STBC does not reduce the ergodic capacity 
of the MIMO channel [13], [22]. For the generator matrix to be unitary, a prerequisite is that the number 
of receive antennas should be atleast equal to the number of transmit antennas, because only then will 
the generator matrix be square. When < nt, only the Alamouti code has been known to achieve 
the ergodic capacity (by saying that an STBC achieves the ergodic capacity, we mean that with the use 
of a suitable outer code in conjuntion with the STBC, capacity can be achieved) of the 2x1 MIMO 
channel. Since it is difficult to make an exact analysis of the ergodic capacity when < nt, we make 
an approximate analysis in the low and high SNR range. The ergodic capacity with the use of a space 
time code is given as follows [22]. 



—tulogdet hn^r H 

2T V nt 



C = —Snlogdet ( hn^r + — HegH^, ) (4) 
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A. Low SNR analysis 

Let HeqH^g = UDU^ be the singular value decomposition of HegH^^. Let D = diag[di ,d2, - ■ ■ , d2Tn^ 
and Heg = [hi,h2, • • • ,h2TnJ- We have, 

1 / SNR r 

2T \ nt 

1 ^ , , SNR^ 

= -^rp^tulogdet hn^T H D 

2T \ nt 

2Trir 



\ i=l 

2Tnr 



nt 



SNR , 

di 



i=l 



f^^H(||He,f) 



i=l 



2ntT 

Since hj = (It ig) H)z;ec(Aj), we have 

\\hif = \\UAif = tr(HAiAf H^). (5) 
The ergodic capacity of an nt x MIMO channel is given as 

C AT p> 

C„,xn. = Siilogdetilrr^ + HH^). (6) 

nt 

In the low SNR scenario, 

Cn,xn. ^ ^^HdlHf). (7) 

nt 

Hence, in the low SNR scenario, if AjAf^ = Vi = 1, 2, • • • , 2Tnr, then, C = Cntxn^, which is 

evident from (O. 
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B. High SNR analysis 

For this purpose, we use the QR decomposition of Hgg. Q and R have the general form obtained by 

Gram — Schmidt process as 

Q = [qi q2 (l3■■■^2Tn^, 



where qj, i = 1, 2, • • • , 2Tnr are column vectors, and 

Ikill (qi,h2) (qi,h3) 
||r2|| (q2,h3) 
R= llrall 



•• (qi,h2Tn, 

•• (q2)h2rn. 



|r2Tn., 



where 



ri = hi, = Ti = hi - Ei=i(qj' hi)qj and q. = i = 2, 3, • • • , 2Tnr. We have, 



C 



Tr\T 



—Sjilogdet hn^T H Hp„H 

2T \ nt 

1 / SNR 

—Snlogdet hn^r + QRR' Q 

21 \ rit 

—Snlogdet hn^T + RR^ 

21 \ nt 

—Syilogdet RR^ 

2T \ nt 



nrlog 



SNR\ 1 



nt 



j + —Snlogdet (RR^) . 



Using the well known identity that the determinant of a triangular matrix is the product of its diagonal 
elements, we have 

, /SNR\ 1 ^ , ,2 
C ^ nrlog { — — ) + —Salog [[ R[t,i) 



\ nt J 2T 



i=l 



2Tnr 



, (SNR\ 1 , 
nrlog I I + —£n ^ logR[ 



1,1) 
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From the definition of the R-matrix, we have, 



Hence, 




Equation ([8]l tells us that at high SNR, the entries of the R-matrix, i.e, (q^, hj) dictate the ergodic 
capacity. If the number of zero entries in the upper block of the R-matrix is larger, then the ergodic 
capacity will be higher. Hence, it is essential that the R-matrix has as many zeros as possible. 

In [5] [Thm. 1], it has been shown that if AjAj^ + = 0„j, then, the i*'* and the j*'^ columns of 

Heg are orthogonal. From the definition of R-matrix, column orthogonality of Hgg dictates the presence 
of zeros. Hence, to design a good STBC when < nt, the equivalent channel matrix should have groups 
of columns orthogonal to one another. We would, of course, like all the columns to be orthogonal, but 
there is a limit to the number, the limit being the maximum number of Hurwitz-Radon matrices for 
nt transmit antennas. Except for the Alamouti code, this number is much lesser than 2Tnr, which is 
the number of weight matrices of a full-rate STBC when < nt. So, evidently, when < nt, higher 
ergodic capacity at high SNR means lower ML-decoding complexity (because of column orthogonality). 
Hence, to construct an STBC with high ergodic capacity, we first construct rate-l STBCs with the lowest 
possible ML-decoding complexity. So far, the known least ML-decoding complexity rate-l codes are the 
rate-l, 4-group decodable codes. But the codes mentioned in literature [15], [16], [17] are not suitable 
for extension to higher number of receive antennas, since their design is obtained by iterative methods. 
Hence, in the next section, we propose a new design methodology to obtain the weight matrices of a 
rate-l, 4-group decodable code by algebraic methods for 2'^ transmit antennas. 

IV. Construction of Rate-1, 4-group decodable codes 

We make use of the following Theorem, presented in [16], to construct rate-l, 4-group decodable codes 
for n = 2" transmit antennas. 
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Al = In 


Aa+i 




A(g-l)fc 1 
g 


A2 


Afc_^2 ~ A2Afc^]^ 

g ^ g 




A(g-l)fc 1 r, = A2A(g-l)fc | -j^ 











Afc 

g 


2k .A/c.Afc_^-j^ 

g g g 




Ak = Afc A(g-i)fc 1 



iABUi i 

Weight matrices of a o-group decodable code 



Theorem 1: [16] An n x n linear dispersion code transmitting k real symbols is (7-group decodable if 
the weight matrices satisfy the following conditions: 

1) A2 = I„,iG{l,2,--- ,|}. 

2) A2 = -L„,jG{2^ + l,m = l,2,... ,<7-l}. 

3) A,A, =A,Aj,i,j G{1,2,--- , |}. 

4) A,Aj=AiA„iG{l,2,--- , |}, j g + 1, m = 1, 2, • • • ,^-1}. 

5) AiAj = -AiAj,i,jG{^^ + l,m = l,2,--- 

6) A^+j = AiA™fc^i,m G {1,2,--- , g - 1}, i G {1, 2, • • • ,g]. 

g g 

Table U illustrates the weight matrices of a (7-group decodable code which satisfy the above conditions. 
The weight matrices in each column belong to the same group. 

In order to obtain a Rate-1, 4-group decodable STBC for 2" transmit antennas, it is sufficient if we 
have 2"+^ matrices satisfying the conditions in Theorem [U To obtain these, we make use of the following 
lemmas. 

Lemma 1: Consider n x n matrices with complex entries. If n = 2*^ and n x n matrices Fj,i = 
1, 2, • • • ,2a anticommute pairwise, then the set of products Fj^Fj^ • • • Fj^ with 1 < ii < • • • < < 2a 
along with I„ forms a basis for the 2^" dimensional space of all n x n matrices over C. 

Proof: Available in [19]. ■ 

Lemma 2: If all the mutually anticommuting nx n matrices Fj, i = 1, 2, • • • , 2a are unitary and anti- 
Hermitian, so that they square to — I^, then the product Fj^Fj^ • • • Fj^ with 1 < zi < • • • < < 2a 
squares to (—1) 2 1„. 



March 16, 2010 



DRAFT 



10 



Proof: (F,,Fi,---F,J(F,,F,,---F,J 

= (-ir-i(F2F,,...F,J(F,,F,3...F,J 

= (-1)-1(-1)-2(f2 p2 . . . F,J(F,3F,, • • • FJ 

= (_1)[(s-1)+(^-2)+-1](f2 p2 ...Yf) 
= (-l)^l™. 

Hence proved. ■ 
Lemma 3: Let Fi,i = 1,2, ••• ,2a be anticommuting, anti-Hermitian, unitary matrices. Let J^i = 
{Fi^,Fi^, • • • ,FiJ and = {Fji,Fj^, • • • ,Fj^} with I < h <■■■< is < 2a and 1 < ji <■■■< jr < 
2a. Let \Qi n 1^2! = P- Then the product matrix Fj^Fj^ • • - Fj^ commutes with Fj^Fj^ • • - Fj^ if exactly 
one of the following is satisfied, and anticommutes otherwise. 

1) r,s and p are all odd. 

2) The product rs is even and p is even (including 0). 
Proof: For Fj^ G ili n 1^2. we note that 

(F,,F,,...F,JF,, = (-l)'^-%,(F,,F,,...FiJ (9) 

and 

(F,,F,,---F,JF,, = (-irF,,(F,,F,,...F,J (10) 

otherwise. Now, 
(F,,,F,,...F,J(F,-F,vF,v) 

= (_1)P(-1)(_1)(-P)^(F^.^F,-, • • -F.-JIF^.F,, • • -F.J 
= (-ir-^(F,,F,vF,v)(F,,F,,...F,J. 

case 1). Since r,s and p are all odd, (— = 1. 

case 2). The product rs is even and p is even (including 0). Hence (— = 1. ■ 
Lemma 4: Let Fj,i = 1,2, ••• ,2a be 2" x 2°' unitary, pairwise anticommuting matrices. Then, the 
product matrix F^^F2^ • • • F2^", Aj G {0, 1}, i = 1, 2, • • • , 2a, with the exception of I20, is traceless. 
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Proof: It is well known that tr(AB) = ir(BA) for any two matrices A and B. Let A and B be two 
invertible, n x n anticommuting matrices. So, 

AB = -BA. 
ABA^ = -B. 

tr{ABA~^) = -tr{B). 
tr(A-^AB) = -tr(B) <^ ir(B) = -tr(B). 

.■.tr(B) = 0. (11) 

Similarly, it can be shown that tr{A) = 0. By applying Lemma [3l it can be seen that any product 
matrix F^^Fg^ • • • Fg^", anticommutes with some other product matrix from the set {F^^Fj^ • • • Fg^", Xi S 
{0, 1}, i = 1, 2, 3, • • • , 2a}. Hence, from the result obtained in ([TT] ). we can say that every product matrix 
F^iFj' • • • except is traceless. ■ 

From Theorem [U to get a rate-1, 4-group decodable STBC, we need 3 pairwise anticommuting, anti- 
Hermitian matrices which commute with a group of 2"-~^ Hermitian, pairwise commuting matrices. Once 
these are identified, the other weight matrices can be easily obtained. From [2], one can obtain 2a pairwise 
anticommuting, anti-Hermitian matrices, presented here for completeness. 

Let 



1 




J 




1 




,P2 = 


,P3 = 




-1 




_i o_ 




-1 



and A®" = A(g)A(g)A---(g)A. 

V ' 

m times 

The 2a anti-Hermitian, pairwise anti-commuting matrices are 

Fi = ±jPf , 

F2 = lf"(8)Pi, 

Fs = lf"(8)P2, 
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For a set 5 = {ai, 02, • • • , a„}, define P(5) as 

P(cS)^{a^a^^...a^,A, G{0,1}}. 

We choose Fi , F2 and F3 to be the three pairwise anticommuting, anti-Hermitian matrices (to be placed 
in the top row along with I„ in Table HI Consider the set 5 = {JF^F^jjFQFj, ■■■ , jF2a-2F2a-i,FiF2F3}, 
the cardinality of which is a — 1. Using Lemma |2] and Lemma [3l one can note that the set consists of 
pairwise commuting matrices which are Hermitian. Moreover, one can note that each of the matrices in the 
set also commutes with Fi, F2 and F3. Hence, P(5), which has cardinality 2""^ is also a set with pairwise 
commuting, Hermitian matrices which also commute with Fi, F2 and F3. The linear independence of P(5) 
over C is easy to see by applying Lemma [T] Hence, we have 3 pairwise anticommuting, anti-Hermitian 
matrices which commute with a group of 2"~^ Hermitian, pairwise commuting matrices. Having obtained 
these, the other weight matrices are obtained from Theorem [T] 

A. An example - n = 8 

To illustrate with an example, we consider the case n = 8. Let Fj, i = 1, 2, • • • ,6 denote the 6 pairwise 
anticommuting, anti-Hermitian matrices. Choose Fi, F2 and F3 to be the three anticommuting matrices. 
Let 

5 = {jF4F5,FiF2F3} 

and 

P(cS) = {Is, jF4F5,FiF2F3, jFiF2F3F4F5}. 

The 16 weight matrices of the rate-1, 4-group decodable code for 8 antennas are as shown below. Each 
column corresponds to the weight matrices in a group. Note that the product of any two matrices in the 
first group is some other matrix in the same group. 
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T 

*8 


17 

1^1 


17 




JF4F5 


iFiF4F5 


JF2F4F5 


JF3F4F5 


F1F2F3 


— F2F3 


F1F3 


— F1F2 


jFiF2F3F4F5 


— jF2F3F4F5 


JF1F3F4F5 


— iFiF2F4F5 



B. Coding gain calculations 

Let A(S,S') = (iet(ASAS^), where AS = S — S',S / S' denotes the codeword difference matrix. 
Let Asi = Si — s'^,i = 1,2, ■ , 2nt, where Si and s- are the real symbols encoding codeword matrices 
S and S', respectively. Hence, 



(2nt 2nt N 

i=l m=l J 



/2nt 2nt \ 
\i=lm=l / 



Note that because of the nature of construction of the weight matrices, 

AiAg = AE^+iA^^^,i, m E {1, 2, 3, 4} p e {1, 2, 3}. 
Further, since the code is 4-group decodable, 

^ ^ AsX, + 2 E AsiAsmAiA^, 

p=0 yi=£|±+l i=^+l m=i+l 

Since all the weight matrices in the first group are Hermitian and pairwise commuting and the product of 
any two such matrices is some other matrix in the same group. It is well known that commuting matrices 
are simultaneously diagonalizable. Hence, 



Ai = ED,E^,iG {2,3,--- ,y} 



where, Dj is a diagonal matrix. Since Aj is Hermitian as well as unitary, the diagonal elements of Dj 
are ±1. In addition, from Lemma|4l there is an equal number of 'I's and '-I's. In fact, because of the 
nature of construction of the matrices Fj, i = 1, 2, • • • , 2a, the product matrices FjFj+i, for even i, and 
the product matrix F1F2F3 are always diagonal (easily seen from the definition of Fj, f = 1, 2, • • • , 2a). 
Hence, all the weight matrices of the first group excluding Ai = I„j are diagonal with the diagonal 
elements being ±1. Since these diagonal matrices also commute with F2 and F3, the diagonal entries are 
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such that for every odd i, if the {i,iY^ entry is 1(-1), then, the (z + l,i + 1)*'* entry is also 1(-1, resp.). 
To summarize, the properties of Aj, i = 2, • • • , ^ are Usted below. 

A, = Af (12) 

A? = In. (13) 

Aj(m, n) = 0,m^n (14) 

Ai{j,j) = ±l,j = 1,2,--- ,nt (15) 

tr(A,) = (16) 

AiijJ) = Mj + 1, j + 1), j = 1,3,5,- ■■ ,nt- I (17) 

AiAj = Ak, i,j,ke {l,2,-- - ,^}. (18) 

In view of these properties, 

g ^ (P + I)nt (.P+l)nt ^ (P + I)it 



where, Dj^ = A^ for some A; G {1, 2, • • • , So, 

Tit d 2 

A(S, S') = det( n E(E ^^^•^^^+^)') 

j=l p=0 1=1 

where, dij = ±1 and dij = 1. In fact, dij = Ai{j,j), i = 1, 2, 3, • • • , ^. Hence, 



min(A(S,S')) = (iet 



Hl^d^jAs, 

j=i \i=i 



From (O, 



2 2 

min(A(S,SO) = det(n(E^ 

i(2j-i)Asi)^). (19) 



We need the minimum determinant to be as high a non-zero number as possible. In this regard, let 

W= =di(2j-i),i,j = 1,2,--- ,y (20) 

and 

yp = [y:it£+i, 2/:i±p+2> ■■■ , y ^tip+i) V = W[s:itp+i, s:iiii+2: • • • > s^^iz+iiV ,P = 0' 1) 2, 3. 
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Lemma 5: W as defined in (1201 ) is a unitary matrix. 

Proof: From ( |20l ). it can be noted that the columns of W are obtained from the diagonal elements 
Aj, i = 1,2, - ■ ■ Each element of a column i of W corresponds to every odd numbered diagonal 
element of Aj. Denote the z*^ column of W by Wj. Applying (fTTl ). (fTSl) and ([T6l) in that order, 

(wi,Wj) = tr(AiAj) 
= tr{Ak) 
= Sii 



where 



Hence, W is unitary. 

Substituting in ( fT9l ). we get 



if i^j 

1 otherwise 



min(A(S,S')) = dei [H; 



So, the minimum determinant is a power of the minimum product distance in nt/2 real dimensions. If 
Yp G Z"^, the product distance can be maximized by premultiplying y^ with a suitable unitary rotation 
matrix U given in [23]. This operation maximizes the minimum determinant and hence the coding gain. 
So, the real symbols of the rate-1, 4-group decodable code are encoded by grouping ^ real symbols 
into 4 groups and each group of symbols taking value from a unitarily rotated vector belonging to Zt", 
the rotation matrix being W^U. The ML-decoding complexity of the code is M~t , where M is the size 
of the complex signal constellation. This is because there are nt/2 real symbols per group to be jointly 
decoded and if we assume a complex constellation of size M, the number of complex symbols to be 
jointly decoded is nt/4. 



V. Extension to higher number of receive antennas 

When rir = 1, a rate-1, 4-group decodable STBC is the best full-rate STBC possible in terms of 
ML-decoding complexity and as a result, ergodic capacity. However, when rir > 1, we need more 
weight matrices to meet the full-rate criterion. Let nt = 2"". We know that if Fj,i = 1,2,- •• ,2a are 
pairwise anticommuting, invertible matrices, then, the set T = {F^^Fg^ • • -Fg^", with Aj G {0, l},i = 
1, 2, • • • , 2a} is a linearly independent set over C. Hence, the set Ai = is linearly independent 
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over M. As a result, the elements of A4 can be used as weight matrices of a full-rate STBC for > 1. 
Keeping in view that the ergodic capacity depends on as many non-diagonal entries of the R-matrix being 
zeros, it is important to choose the weight matrices prudently. The idea is that given a full-rate STBC 
for rir — 1 receive antennas, obtain the additional weight matrices of a full-rate STBC for Ur receive 
antennas by using the weight matrices of a rate-1, 4-group decodable STBC such that the after addition 
of the new weight matrices, the set of weight matrices is linearly independent over M. This is achieved 
as follows. 

1) Obtain a rate-1, 4-group decodable STBC by using the construction method detailed in Section 
HVl Due to the nature of construction, the product of any two weight matrices is always some 
other weight matrix of the code, up to negation. Denote the set of weight matrices by Qi. The 
ML-decoding complexity of the code is M^. 

2) From the set T, choose a matrix that does not belong to Qi and multiply it with the elements of 
Qi to obtain a new set of weight matrices, denoted by Q2- Clearly, the two sets will not have any 
matrix in common. The weight matrices of Q2 form a new, rate-1, 4-group decodable STBC. This 
is because the ML-decoding complexity does not change by multiplying the weight matrices of a 
code with a unitary weight matrix. In this case, we have multiplied the elements of Qi with an 
element of T, which is a unitary matrix. Now, Qi [j Q2 is the set of weight matrices of a rate-2 
code with an ML-decoding complexity of M^'.M^ = M^. This is achieved by decoding the 
last rit symbols with a complexity of M"* and then conditionally decoding the first rit symbols 
using the 4-group decodability property. 

3) For increasing number of n^, repeat as in the second step, obtaining new rate-1, 4-group decodable 
codes and then appending their weight matrices to obtain a new, rate-n,. code with an ML-decoding 
complexity of M""^"'-"!). 

4) When all the elements of T have been exhausted (this occurs when = nt/2), the remaining 
matrices up to a rate of nt symbols per channel use can be obtained from jF. Note from Lemma 
[T] that this does not spoil the linear independence over M of the weight matrices. 

The R-matrix of the STBC for rir receive antennas has the following structure, irrespective of the 
channel realization. 



D 



X 



X 



R = 



02„, 



D 



X 



(21) 



02nt 02nt 



D 



Mai-ch 16, 2010 



DRAFT 



17 



where X G ]^2ntx2nt ^ random non-sparse matrix whose entries depend on the channel coefficients 
and D = I4 (g) V, with V G M^^t" being an upper triangular matrix. As a result of the structure of 
D, the R-matrix has a large number of zeros in the upper block, and hence, compared to other existing 
codes, the proposed codes are expected to have higher ergodic capacity (for n,. < n^) and lower average 
ML-decoding complexity. 

A. The Silver code as a special case of rit = 2 

The silver code, which is well known for being a low complexity, full-rate, full-diversity STBC for 
rir > 2, transmits 2 complex symbols per channel use. Its first four weight matrices are that of the 
Alamouti code, which is a rate-1, 4-group decodable STBC for 2 transmit antennas. The Silver code's 
next 4 weight matrices are obtained by multiplying the first four weight matrices with j. However, to 
make the code a full-ranked one, the last four symbols take values from a different constellation, which 
is obtained by unitarily rotating the symbol vector in Z[j]^. The Silver code compares very well with 
the well known Golden code in error performance, while offering lower ML-decoding complexity of 
for square-QAM only. 

VI. Simulation results 

In all the simulation scenarios in this section, we consider the Rayleigh block fading MIMO channel. 
We consider 8 transmit antennas. To construct a rate-2 code for 8 transmit antennas, we first construct 
a rate-1, 4-group decodable STBC as described in Section UVl and denote the set of obtained weight 
matrices by Qi. Next we multiply the weight matrices of Qi by F4 to obtain a new set of weight matrices 
which is denoted by ^2- The weight matrices of the new STBC with rate-2 are obtained from Qi [j Q2- 
A rate-3 code for 3 receive antennas can be obtained by multiplying the matrices of Qi with Fg and 
appending the resulting weight matrices to the set Qi |J Q2- The rival code is the punctured perfect code 
for 8 transmit antennas [14]. The ergodic capacity plots of the two codes are shown in Fig.[T] As expected, 
our code achieves higher ergodic capacity, although lower than that of the corresponding MIMO channel. 
It must however, be noted that both codes help to achieve the same ergodic capacity as that of the MIMO 
channel for > 8 because the generator matrix is unitary in that case. 

Fig. |2] shows the codeword error performance of our code for 8 x 2 system and the punctured perfect 
code using 4-QAM. The performance is more or less the same. Unlike the Perfect code, our code will 
not have full-diversity if the design is made as explained in Section |IVl However, we have multiplied 
the weight matrices of Q2 with the scalar in order to enhance performance. Fig. [3] shows the symbol 
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SNRindb 

Fig. 1. Ergodic capacity Vs SNR for codes for 8 x 2 and 8x3 systems 

eiTor performances of the two codes. Our code has a better performance at a higher SNR. The reason 
for this is that the number of symbol errors per codeword error for the punctured Perfect code is more 
than that of our code. Hence, even though the CER is the same, the SER is different. Our code appears 
to have full diversity, but we have not been able to prove it. The most important aspect of our code is 
that it has an ML-decoding complexity of M^^, while that of the comparable punctured Perfect code is 

VII. Discussion 

In this paper, we proposed a scheme to obtain a full-rate STBC for 2" transmit antennas and any number 
of receive antennas with reduced ML-decoding complexity. The STBCs thus obtained have higher ergodic 
capacity at high SNR than existing STBCs for the case < n^. We have, however, not been able to 
provide a scheme to obtain full-diversity codes from these designs. Also it is to be seen if the proposed 
codes are better suited than existing codes for sub-optimal decoding techniques like lattice reduction 
aided detection, owing to the fact that more number of symbols are disentangled from one another than 
in the case of known codes. These are some of the directions for future research. 
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